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Abstract 

Background: The biomechanical property of the periodontal ligament (PDL) is 
important in orthodontics and prosthodontics. The objective of this study was to 
evaluate the feasibility of measuring the biomechanical behavior of the periodontal 
ligament using micro-computed tomography (micro-CT). 

Methods: A custom-made apparatus measured the force and displacement of a 
porcine PDL specimen within the micro-CT environment. Synchronized computed 
tomography (CT) images were used to obtain the deformation and displacement of 
the entire specimen and to reconstruct the three-dimensional mesh model. To 
match the experimental results, finite element analysis was then applied to simulate 
the biomechanical response of the PDL. The mechanical model of the PDL was 
assumed as the hyperelastic material in this study. 

Results: The volume variations of the tooth and the alveolar bone were less than 
1%, which implies that tooth displacement was caused mostly by displacement of 
the PDL. Only translational displacement was observed with each load step because 
the transformation matrix acquired from the CT image registration was identical. The 
force-displacement curve revealed the nonlinear behavior of the PDL. There was a 
high correlation between the experimental displacement results and the simulation 
displacement results. The numerical results (based on the assumption that the PDL is 
the hyperelastic material) showed good agreement with the experimental results. 

Conclusions: Nondestructive measurements by micro-CT obtained the biomechanical 
behavior of the PDL. Using the hyperelastic characteristic as the constitutive model can 
properly predict the force-displacement relation of the PDL after loading. This study 
provided a feasible approach for measuring the biomechanical behavior of the PDL for 
further dental application. 
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Background 

The periodontal ligament (PDL) has a determinative role in dental biomechanics. The 
PDL connects the tooth and alveolar bone, absorbs occlusal impact, and forms and re- 
sorbs the alveolar bone for tooth movement. Hence, the biomechanical property of the 
PDL is important in orthodontics and prosthodontics. Orthodontic treatment is a re- 
iterative and time-consuming process during which a tooth gradually moves to its 
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prospective position. The physiological mechanism of tooth movement is primarily be- 
cause of the response of the PDL while an orthodontic force is applied [1]. 

The PDL is clearly important, but its biomechanical behavior remains unclear. Toms 
et al. [2] performed a shear test by using a transverse thin specimen of the premolar 
tooth. They revealed the anisotropic behavior of the PDL. Dorow et al. [3] designed an 
apparatus to test the response of the porcine anterior tooth under different strain rates. 
Their results showed that the nonlinear hysteresis behavior varied with the strain rate. 
Natali et al. [4] designed a laser-optical system to measure the minipig premolar tooth. 
They obtained the nonlinear force-displacement relation of the PDL. In 2009, Tohill 
et al. [5] acquired the relaxation characteristic of the premolar porcine tooth. However, 
these studies only demonstrated the simplex mechanical feature of the PDL. 

To date, the PDL cannot be harvested by a nondestructive approach for evaluating 
the stress-strain relation. Analytical models have consequently become a popular ap- 
proach for estimating the biomechanical behavior of the PDL. There are three types of 
governing equation for these models: (1) linear elastic model, (2) viscoelastic model, 
and (3) hyperelastic model. The linear elastic behavior of the PDL has been widely in- 
troduced in previous literature reports, but the Youngs modulus values range from 
0.07 MPa to 1750 MPa [6]. However, some studies indicate that the elastic property of 
the PDL is nonlinear [7-9]. Previous analytical models have assumed viscoelastic prop- 
erties for the PDL such as creep, stress relaxation, and hysteresis [10-13]. However, 
most PDL researchers prefer the linear elastic behavior because of its simplicity. Be- 
cause the force-displacement relation of the PDL is insufficient, an energy-based equa- 
tion was developed and named the "hyperelastic model". To predict tooth mobility or 
movement, many researchers have tried to establish a hyperelastic model of the PDL 
[4,13-15]. Choosing hyperelastic models and parameters is arduous and may directly in- 
fluence the simulation results. The dilemma of predicting PDL behavior has limited the 
advancement of dental biomechanics. 

Micro- computed tomography has recently been widely applied for dental research in 
the microscale such as dental morphology [16], dental biology [17], dental biomechan- 
ics [18,19], and dental material science [20]. This method can detect the precise 
morphology and can inspect the details of the microstructure or microvariations (e.g., 
movement or deformation) in biological tissues. Therefore, we aimed to measure the 
biomechanical behavior of the PDL using micro-CT. 

Methods 

Figure 1 shows our research protocol. A porcine premolar tooth with the mandible was 
harvested for biomechanical testing. A device was designed to synchronously measure 
the reaction force under controlled displacement in the micro-CT environment. The 
sectional micro-CT images were simultaneously digitized for three-dimensional (3D) 
FE model reconstruction. The force-displacement diagrams of the FE and the experi- 
mental results were compared to determine the material parameters of the PDL. 

Custom device for installation into the micro-CT scanner 

For simultaneously obtaining the reaction force and the displacement of the entire pre- 
molar tooth, a device was developed for the micro-CT environment. The device was 
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Determination of PDL Force-Displacement Diagram 

Material Parameters 

Figure 1 The research protocol of the study. 



able to control the displacement and to measure the reaction force. As shown in 
Figure 2, the main components of the device were (1) a micrometer (Model No. 103- 
137, Mitutoyo Corporation., Japan) with a sensitivity of 0.01 mm to control the dis- 
placement of the premolar tooth in our experiment; (2) a universal joint (JA20-8-125 
floating joint; SMC Corporation, IN, USA) to ensure that the direction of the reaction 
force would remain axially loaded; (3) a load cell (MDB-50; Transducer Techniques, 
CA, USA) to record the magnitude of the axial force; (4) a linear motion rolling guide 
(IKO LWL9 B; Nippon Thompson Co., Ltd., Japan) in both directions: first, to maintain 
the same direction between the load cell and the micrometer and, second, to drive the 
sample jig to apply tension or a compression force to the specimen; and (5) the sample 
jig to fix the specimen and apply force. The entire device was designed based on the 
folder size of micro-CT. It could record the force applied to the tooth. 



Experiment 

A porcine premolar tooth was harvested, but the surrounding gingiva and mucosa were 
eliminated. The complete specimen included the premolar tooth, PDL, and alveolar 
bone. The crown and alveolar bone were both perforated with a 3-mm hole to fix them 
into the sample jig. The premolar tooth was then loaded by micrometer with constant 
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displacements of 0.17 mm, 0.34 mm, 0.60 mm, 0.84 mm, and 1.28 mm while the reac- 
tion force was recorded in the end of the micro-CT scan procedure. Sectional images 
of the entire specimen were scanned by micro-CT (SkyScan 1076; Bruker-MicroCT, 
Belgium) for further image processing and for the reconstruction of the finite element 
model. The parameters of the micro-CT were 70 kV, 100 (lA, 35- (im resolution, and 
316-ms exposure time. 

The real displacement of the premolar tooth was calculated from the preloaded and 
postloaded stereolithography (STL) model that was reconstructed from micro-CT im- 
ages by medical image software (Mimics vlS.Ol; Materialise, Belgium). The preloaded 
alveolar bone and tooth models were respectively reconstructed as the references. The 
perforated holes in the alveolar bone and tooth were selected as the registration points. 
By superimposing the postloaded alveolar bone model, the displacement of the tooth 
could be obtained (Figure 3). During image processing, the alveolar bone was assumed 
to be a rigid body (i.e., the alveolar bone had no deformation or displacement). A built- 
in function (i.e., STL registration) was adopted to compare the volumes of the alveolar 
bone and tooth to ensure that variation was mostly caused by the displacement of the 
tooth, rather than by bone deformation. In addition, a 4 x 4 transformation matrix 
could be obtained as the real displacement of the premolar tooth. 



Finite element analysis 

After building a solid model of the specimen from the STL file, the 3D finite element 
model (which included the premolar tooth, pulp, PDL, and alveolar bone) was 




Figure 3 The STL registration results of the alveolar bone and tooth (applied displacement from 
left to right: 0.17 mm, 0.34 mm, 0.60 mm, 0.84 mm and 1.28 mm). (A) Colored contour represented 
that the postloaded alveolar bone was superimposed to the preloaded bone. (B) Colored contour 
represented that the postloaded tooth was superimposed to the preloaded tooth. 

V ) 
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developed in the FE package (ANSYS; ANSYS, Inc., Canonsburg, PA, USA) (Figure 4). 
A 10-node tetrahedral element (Solid 187) was used in this study. The material proper- 
ties of the premolar tooth (E = 22000 MPa; v = 0.3) and alveolar bone (E = 1200 MPa; v 
= 0.3) were defined as isotropic, linear elastic [14]. The material property of the PDL 
was assumed to be hyperelastic (i.e., a 3 parameters Mooney-Rivlin model). The fol- 
lowing formula was used: 

W = Cio(/i-3) + Coi(/2-3) + Cn(/i-3)(/2-3) +^(/-l)' 

in which Qi, and Cn are the material constants and d is the material incom- 
pressibility parameter. Ii and [2 are the first and second deviatoric strain invariant. The 
initial estimation of these material constants (Q^ = 0.04 MPa; Coi = 0.02 MPa; Cn- 
0.04 MPa; and <i = 0.02), to which previous researchers have referred [11], was estab- 
lished to obtain a proper combination for fitting the experiment results. 

The boundary and loading conditions were followed by the experiments. The surface 
nodes of the perforated hole in the alveolar bone were fixed in every direction and the 
loadings, according to the experimental results, were applied to the surface nodes of 
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the perforated hole in the tooth. The all interfaces between tissues were continuous to 
ensure the loading could be properly transmitted. 

Results 

Table 1 showed the alveolar bone and tooth volumes at each load step. Volume change 
variations, which were reconstructed and computed from the serial CT images by STL 
registration, were all less than 1%. This result implied that the alveolar bone and tooth 
remained mostly undeformed during loading (i.e., the displacement was primarily 
caused by the PDL). In addition, a 4 x 4 transformation matrix (representing the trans- 
lational and rotational information of the tooth) could be acquired from the STL regis- 
tration. The results of the transformation matrix at each load step approximated the 
identity matrix which the main diagonal elements are Is and all the remaining ele- 
ments are 0 s. This result implied that a rotational effect could be ignored, whereas the 
translational term in the occlusal-gingival direction of the transformation matrix could 
be treated as the real displacement of the tooth. The actual displacement of the tooth 
could be obtained from the transformation matrix. These results demonstrated that the 
PDL nearly possesses the bilinear characteristic of stiffness, which increases with the 
applied force. 

The material constants Cio, Cqi, Cjj^ and d were obtained after trial-and-error FE 
analysis. Through the trial-and-error process, we found that the constants Cji and d 
had a greater effect on the biomechanical properties of the PDL. At the end of the ana- 
lysis, all material constants were determined as follows: €20= 1 x 10"^ MPa; C^i = 1 x 
10"^ MPa; Cjj = 0.1 MPa; and d = 4,n gure 5 demonstrates the linear regression of the 
experimental displacement and the FE estimated displacement. The coefficient of deter- 
mination was 0.98. 

Discussion 

In this study, micro CT was demonstrated to be a useful tool for measuring the bio- 
mechanical behavior of the periodontal ligament. This approach can be used to deter- 
mine the response of the PDL according to the entire displacement with external 
stimuli applied to the tooth. Many studies have investigated the biomechanical behavior 
of the PDL by using specimens from animals such as rats, rabbits, and pigs. [4,7,21,22]. 
Animal specimens are used because intact specimens of the tooth and adjacent alveolar 
bone from human cadavers are difficult to harvest. This study preserved the intact por- 
cine PDL and adjacent bone for the investigation of PDL biomechanics. Furthermore, 

Table 1 The variations in the volume of the alveolar bone and the tooth at each load 
step (relative to zero force and displacement) 



Applied Displacement Alveolar Bone Volume (mm^) Error (%) Tooth Volume (mm^) Error (%) 



0 mm 


4117.65 




322.14 




0.17 mm 


4085.32 


0.7 


323.93 


0.6 


0.34 mm 


4080.62 


0.8 


324.01 


0.6 


0.60 mm 


4088.57 


0.7 


323.24 


0.3 


0.86 mm 


4080.59 


0.9 


323.33 


0.4 


1.28 mm 


4106.22 


0.3 


322.53 


0.1 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 

FE Displacement (mm) 

Figure 5 Linear regression of the experimental (i.e., actual) displacement and the finite element (FE) 
estimated displacement. 



using the custom-made apparatus in the micro-CT environment made it possible to 
simultaneously record the displacement and the corresponding force of the tooth and 
the PDL. This setup may create an in vitro environment that is biomechanically com- 
patible with the loading experienced by the tooth and PDL during occlusion. 

This study applied a hyperelastic material model to simulate the mechanical behavior 
of the PDL. The displacement results between the in vitro experiment and the finite 
element analysis showed a high correlation (Figure 5), which implied the force- 
displacement results from FE analysis were similar to that from the experiment. This 
study provided a feasible approach to evaluate the holistic material model of the PDL. 
The difference between the estimated displacement and the experimental displacement 
was less than 15%. The force-displacement results also demonstrated that the PDL ap- 
proximately obeyed hyperelastic behavior (Figure 6). However, the coefficients used in 
this study were incongruous to those reported in previous studies [11]. This incongruity 
may be because the experimental and finite element models used in this study were 
based on the intact PDL of the premolar tooth. Besides, only one tooth was adopted in 
this study for comparison. Nevertheless, micro CT still provided a feasible modality to 
measure the biomechanical behavior of the PDL. 

Although it is easier to apply a simple linear elastic model, this assumption was not 
suitable for the soft tissue, which was well-known as a nonlinear material. The general 
behavior of the soft tissue was divided into three regions: toe, linear and failure regions. 
When the loading of the soft tissue increased initially, the stiffness of the toe region 
was low; the linear region was identified after the toe region and had a higher stiffness 
than the toe region; then the soft tissue failed eventually. It is not easy to clarify the 
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transition zone between each region. Furthermore, the viscoelastic model was also not 
suitable in this study because the experimental results did not reveal the time- 
dependent characteristic of the PDL. Therefore, the hyperelastic model was more favor- 
able to fit the nonlinear force-displacement behavior of the PDL. 

The results of previous studies represented the local characteristic of the PDL instead 
of the whole tooth structure; hence, the biomechanical behavior would greatly differ 
among these studies. Therefore, finding a suitable constitutive model of the PDL is an 
important issue in dental biomechanics, especially in the development of prosthodon- 
tics and orthodontics. However, the parameters reported in this study resulted in 
greater error when a small force was applied. In addition, each tooth might have its 
own constitutive equation according to the morphology; therefore, the parameters may 
still need further validation. 

From the sectional CT images, we observed that the transversal thickness of the PDL dif- 
fered, based on the location of the tooth. The transversal thickness of the PDL was greater 
around the alveolar crest and the apical region than in other regions. Previous studies have 
also found similar results [23,24]. In addition, most specimens were sliced transversely into 
a thin plate for performing a shear test [24], which could only represent the local biomech- 
anical behavior of the PDL. Therefore, this study demonstrated the global biomechanical 
behavior of the PDL and may have greater differences in comparison to other studies. Fur- 
thermore, because multiple roots of a tooth could disturb the biomechanical performance 
of the PDL during occlusal pressure, this study applied a tension force to highlight the inde- 
pendent effect of the PDL. The resistance force increased dramatically with the 
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displacement after passing the toe region, based on the experimental results (Figure 6). This 
would explain the protective mechanism that prevents the tooth from excessive extrusion. 

Testing in the micro-CT environment has several advantages compared to the previous 
experiments. First, it is a nondestructive modality for performing biomechanical testing. 
Further biological or physiological analysis can be performed because the PDL specimen re- 
mains intact. Second, micro or local observations can be accomplished because micro-CT 
has a higher resolution than traditional modalities such as CT, magnetic resonance imaging, 
or ultrasound images. Analysis of the intact PDL or observation of internal morphology is 
available by 3D model reconstruction because no permanent damage occurs within the 
physiological loading range. 

Bone remodeling is an important process in orthodontic biomechanics. Studies report 
a direct correlation of stress-strain fields in the PDL with alveolar bone resorption 
[25,26]. While PDL can be a direct effector, the response of bone to mechanical load 
could be a delayed effect that is necessary for tooth translation. Instead of a time- 
dependent (i.e., viscoelastic) characteristic of the PDL, our results showed good agree- 
ment between the numerical simulation and the experiment. This implied that the 
hyperelastic constitutive model is practical for investigating the biomechanical response 
of the PDL. The mechanobiological interaction of the PDL and alveolar bone neverthe- 
less remains indeterminable. 

Our results demonstrated that the PDL is suitable for simulating hyperelastic mater- 
ial; however, some limitations still need to be depicted. Because of the scattering effect 
of the micro-CT, the material selected for the custom-made apparatus should be taken 
in account to eliminate its influence in image registration. The experimental results 
and the hyperelastic model can well interpret the steady state of the PDL after loading 
because the entire scanning process requires more than 30 minutes. More samples 
under various loading modes are required for a complete evaluation of the biomechan- 
ical property of the PDL. The detailed tooth structures were not modelled in this study, 
such as enamel, dentine and pulp. The Young s modulus and stiffness of the anatomical 
layer of enamel and dentine was much greater than those of the PDL. In addition, the 
volume of pulp took a very small part of the entire tooth. Furthermore, from the results 
of micro-CT reconstruction, the tooth and bone was approximately undeformed. 
Therefore, the lack of these layers in FE model would not affect the results much. This 
study excluded the behavior of the PDL at a high strain rate. For this study, only one 
specimen of a porcine mandibular premolar tooth was harvested. However, finding a 
suitable constant set of the hyperelastic model was a time-consuming task. Besides, the 
individual difference was existed in every sample harvested from the porcine. Although 
more specimens may have enhanced the reliability of the experimental and numerical 
results, this study successfully established a feasible approach to predict the force- 
displacement relation of the PDL in a holistic observation. 

Conclusion 

Micro-computed tomography, which performs nondestructive measurements, is cap- 
able of obtaining the biomechanical behavior of the PDL. Using the hyperelastic charac- 
teristic as the constitutive model may properly predict the force-displacement relation 
of the PDL after loading. 
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